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FORWORD 

This  report  summarizes  the  methods  and  key  findings  of  a 
fluid  modelling  study  performed  in  the  BBN  Atmospheric  Boundary 
Layer  wind  tunnel  to  simulate  the  fluid  dynamic  dispersion  of 
pollutants  released  from  two  large  ventilation  stacks  associated 
with  the  proposed  Copley  Place  project.   The  facility  and  model 
used  were  identical  to  that  described  in  Appendix  A  of  a  companion 
report  entitled  "Pedestrian  Wind  Environment  at  Copley  Place, 
Boston,  Massachusetts"  (BBN  Report  Number  4415) . 

The  simulation  of  such  dispersion  processes  involves  a 
reasonably  elaborate  documentation  of  flow  phenomena  and  a  care- 
ful monitoring  of  tracer  gas  concentrations  released  into  the 
system.   The  methods  used  for  performing  this  work  are  briefly 
described  in  the  report,  and  are  consistent  with  those  recommended 
in  the  EPA  draft  guideline  for  fluid  dynamic  modelling*,  which 
is  available  for  public  inspection. 

The  results  of  this  study  are  presented  in  the  form  of  dilu- 
tion factors  which  related  to  concentration  of  pollutants  at  any 
point  outside  of  the  stack  to  that  in  the  stack.   These  dilution 
factors  can  then  be  applied  to  the  full  scale  system  for  any 
pollutant   emission  rate  from  within  the  tunnel  or  garage  areas. 
The  incremental  increase  in  the  air  pollution  may  then  be  added 


*Snyder,  W.H.  "Guidelines  for  Fluid  Modelling  of  Atmospheric 
Dispersion"  EPA-450/4-79-016  June  1979  (Draft) 
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to  the  existing  or  predicted  background  levels  from  other  sources 
such  as  street  vehicular  traffic,  and  the  hourly  or  eight-hour 
CO  concentrations  can  be  predicted  and  compared  with  standards. 
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I.   INTRODUCTION 

The  Copley  Place  project  will  cover  the  Mass.  Turnpike  east 
and  west  lanes,  the  B  &  A  railroad  line,  and  several  turnpike 
ramps.   In  addition,  the  project  will  contain  two  parking  garages 
which  must  be  ventilated  to  maintain  favorable  air  quality  within 
the  garages.   Two  ventilation  systems  are  to  be  operated  separately 
using  fans  to  push  air  through  ventilation  stacks,  which  are 
located  along  the  south  side  of  Stuart  Street  (Figure  1) .   The 
ventilation  stack  for  turnpike  ramps  B  and  D,  and  the  Marriott 
and  Central  garages,  is  located  in  the  southeast  corner  of  the 
retail  area  approximately  at  roof  level  (elevation  185)  ,  also 
shown  in  Figure  1.   The  design  flow  rate  for  these  two  stacks  at 
full  capacity  are  1.23  million  cubic  feet  per  minute  (cfm)  for 
turnpike  ventilation  stack  (Stack  I)  and  640,000  cfm  for  the 
ramp  and  garage  ventilation  stack  (Stack  II) . 

Under  certain  wind  directions  and  speeds,  the  location  of 
these  ventilation  stacks  within  the  project  elements  gives  rise 
to  aerodynamic  phenomena  (called  downwash)  which  cause  the  fumes 
and  pollutants  being  ventilated  through  these  stacks  to  be  brought 
down  to  street  level  in  areas  near  high-rise  and  other  buildings. 
Likewise,  the  "canyoning"  effect  of  the  various  streets  can  reduce 
the  rate  of  dispersion  of  air  pollutants  from  these  sources  below 
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Figure  1.   Location  of  Ventilation  Stacks  and  Receptors  within 
Copley  Place  Site. 
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what  it  would  be  in  the  absence  of  large  street  "canyons".   The 
ability  of  present  day  air  pollution  dispersion  computational 
models  to  accurately  predict  the  dispersion  phenomena  associated 
with  sources  released  near  high-rise  elements  or  in  complex 
street  "canyons"  is  very  limited.   The  consequence  of  the  limita- 
tions of  computational  models  may  either  be  an  over-estimate  or 
an  under-estimate  of  the  pollution  levels  at  various  points  near 
the  Copley  Place  Project  or  similar  projects.   The  consequence 
of  an  overestimate  would  lead  to  design  of  an  unnecessarily 
elaborate  and  expensive  ventilation  system,  while  an  underestimate 
of  the  levels  could  lead  to  localized  areas  of  unacceptably  high 
levels  of  carbon  monoxide  or  other  fumes. 

Wind  tunnel  studies  of  these  dispersion  processes  have  led 
to  improvements  in  the  accuracy  of  prediction  of  air  pollution 
impacts  in  complex  urban  environments.   The  next  section  of  this 
report  briefly  summarizes  the  techniques  which  are  used  to  conduct 
such  a  modelling  study;  Section  3  summarizes  key  results. 
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2.   SUI4I1ARY  OF  METHODOLOGY 

The  modelling  of  dispersion  from  point  sources  in  complex 
environments  involves: 

(11   Proper  simulation  of  the  geometry  of  the  area, 

(2)   The  detailed  correct  simulation  of  the  mean  velocity 
profile  and  atmospheric  turbulence  characteristics 
incident  upon  the  site  and  within  the  stack  emitting 
the  pollutant  of  interest.   This   simulation  of  turbulence 
and  flow  parameters  for  both  the  atmosphere  and  the 
source  are  crucial  to  the  proper  achievement  of  meaningful 
results  within  the  wind  tunnel; 

C3)   A  simulation  of  both  buoyancy  and  momentiom  effects  within 
a  stack.   If  the  stack  gas  is"cold"and  consists  largely 
of  air  at  ambient  temperatures  with  small  amounts  of 
pollutants  or  other  concentrations,  then  the  buoyancy 
effects  of  the  plume  emanating  from  the  stack  are 
negligible  and  the  momentum  of  the  jet  dominates  the 
mixing  of  the  pollutants  with  ambient  air  and  thus  the 
dispersion  process. 

(4)   Measurement  of  ground  level  locations  using  a  system  which 
can  resolve  pollution  at  small  points  without  disturbing 
the  flow  field  in  the  neighborhood. 
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The  methods  used  in  the  BEN  atmospheric  boundary  layer  wind 
tunnel  for  fluid  dispersion  modelling  are  consistent  with  those 
recommended  by  EPA  in  a  recently  proposed  guideline  document.  * 
For  the  present  study  it  was  determined  that  the  so-called  Froude 
number  (ifatio  of  inertial  forces  to  buoyant  forces  in  a  plume) 
would  be  simulated  without  artificially  lightening  the  gas  in  the 

stack  through  the  use  of  helium  or  other  lighter-than-air  gases. 
The  dispersion  processes  are  governed  by  turbulent  mixing  of  the 
plume,  which  is  due  to  the  side  shear  layers  of  the  plume  and 
atmospheric  turbulence.   This  process  is  quite  dependent  upon  a 
dimensionless  ratio  of  fluid  inertial  forces  to  those  due 
to  the  internal  viscosity  of  the  fluid  ;  this  parameter  is  known 
as  the  Reynolds  Number.   It  is  not  always  possible  in  model 
scale  to  duplicate  the  full  scale  Reynolds  number  to  a  degree 
sufficient  to  achieve  identical  mixing  processes  without  artifi- 
cially disturbing  the  flow.   Thus,  the  flow  in  the  jet  is  made 
to  behave  like  a  full-scale  flow,  i.e.  one  at  full  scale  Reynolds 
Number, by  artificially  introducing  turbulence  in  the  jet  through 
the  use  of  screens  and  other  flow  disturbing  devices.   The 
consequence  of  a  low  Reynolds  number  is  less  turbulent  mixing 
than  that  at  full  scale^  and  thus  higher  local  concentrations 

The. methodology  used  to  acquire  data  in  the  wind  tunnel  is 
now  described.   The  system  used  consists  of  a  flow  mixing  system 


*  See  earlier  footnote 
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which  supplies  air  to  the  exhause  stack  (or  stacks) ,  which  is 
mixed  with  a  tracer  gas,  methane  (CH^j ,  in  relatively  low  concen- 
trations, such  that  the  mass  flow  in  the  stack  is  dominated  by 
the  clean  air  system.   The  stack  exit   is  calibrated  to  ensure 
uniformity  of  velocity  profile  and  methane  concentrations. 

The  wind  tunnel  is  then  operated  at  a  variety  of  speeds 
to  enable  examination  of  the  full  range  of  dispersion  parameters. 
Receptor  nozzles  consisting  of  small  tubes  anchored  at  various 
positions  on  the  model  are  placed  at  so-called  key  receptor 
locations  at  various  points  around  the  site  and  the  neighborhood 
adjacent  to  the  project.   These  nozzles  may  be  either  sampled 
continuously  to  obtain  a  time  history  of  the  concentration  at  a 
given  point,  or  by  a  simultaneous  sampling  system  which  acquires 
and  holds  samples  from  each  of  up  to  80  nozzles  simultaneously. 
Both  capabilities  of  the  nozzle  sampling  system  may  be  useful 
in  a  project  such  as  this.   However,  to  enable  the  experimenters 
to  quickly  locate  wind  azimuths  and  receptors  which  may  be  critical 
to  the  environmental  impact  of  the  ventilation  system,  the 
receptors  are  operated  in  a  continuous  mode  one  at  a  time.   The 
output  of  the  sampling  system  is  continuously  plotted  in  time 
while  the  turntable   holding  the  model  is  slowly  rotated  through  360 
degrees  of  wind  direction.   A  sample  plot  which  results  from  this 
process  is  shown  in  Figure  2.   The  traverse  rate  on  this  plot 
is  between  one-half  degree  and  one  degree  of  wind  azimuth. 
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per  second.   This  causes  some  rapid  fluctuation  in  the  instant- 
aneous concentration  received  at  the  receptor;  however,  these 
fluctuations  are  related  to  aerodynamic  phenomena  which  are 
slowly  changing  as  the  model  rotates  and  are  quite  typical  of 
the  actual  situation  to  be  expected  at  full  scale. 

Upon  identifying  the  wind  azimuths  at  which  peak  concentrations 
occur,  the  model  is  returned  to  those  wind  directions,  and  the 
nozzle  is  sampled  for  a  longer  period  of  time  to  measure  a  long- 
term  average  ,  and  standard  deviation  from  the  average,  for  the 
flow  conditions  in  interest.   Also  at  this  point,  variations  in 
wind  velocity  and  stack  effluent  velocity  are  made  to  determine 
the  parametric  behavior  of  the  concentration  reaching  a  parti- 
cular receptor.   A  sample  result  of  those  parametric  variations 
is  shown  in  Figure  3,  for  a  receptor  located  along  Stuart  Street 
facing  a  low  rise  element  of  Copley  Place  Project,  on  the  (north) 
Western  International  side  of  Stuart  Street.   The  source  of 
pollutant  in  this  case  is  Stack  II,  operated  at  full  load,  1/2 
load,  and  1/4  load,  all  of  which  are  possible  options  in  the  opera- 
tion of  the  actual  system.   It  can  be  seen  that  at  very  low  wind 
speeds,  the  plume  apparently  escapes  any  influence  of  the  buildings 
and  mixes  with  the  atmosphere  at  very  high  altitudes,  never 
reaching  the  ground.   As  the  wind  speed  increases,  the  plume  is 
bent  over  and  impinges  on  the  tops  of  buildings  and,  subsequently 
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becomes  entrapped  in  the  downwash  between  the  Western  International 
Hotel  and  the  low  rise  along  Stuart  Street;  this  entrapment  leads 
to  a  peak  reading  on  the  Stuart  Street  side  which  in  this  case  is 
around  0.2  parts  per  million  of  carbon  monoxide,  when  the  results 
are  scaled  to  full  scale  conditions*.   It  can  be  seen  that  when 
the  system  is  operated  at  full  flow, the  magnitude  of  this  peak 
concentration  can  be  reduced.   It  is  also  seen  that  as  the 
gradient  level  wind  speeds  increase,  the  atmospheric  mixing 
processes  tend  to  reduce  the  concentration  of  pollutant  reaching 
the  street.   (The  mean  annual  wind  speed  for  Boston  is  between 
28  and  30  fps  at  gradient  height;   thus,  for  average  day  conditions, 
the  stack  plumes  will  be  well-mixed  by  the  atmospheric  winds.) 
Similar  parametric  curves  are  shown  for  receptor  E,  which  is  at 
the  corner  of  the  southwest  corridor  deck  and  Harcourt  Streets, 
for  the  Stack  II  operation  (Figure  4) ,  and  Receptor  A  when  the 
larger  Stack  I  is  operated  (Figure  5) .   Receptor  A  receives  the 
largest  concentration  of  ground  level  pollutants  form  the  venti- 
lation system  of  any  examined  in  the  entire  system. 


*See  Appendix  A  for  scaling  relationship 
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3.   SUMMARY  OF  RESULTS 

The  results  of  these  studies  can  be  summarized  in  the  form 
of  dilution  factors  for  the  various  combinations  of  stack  flow  , 
wind  speeds,  and  wind  directions  at  each  receptor  and  can  also 
be  translated  to  full  scale  CO  concentration  increases  (above 
background  level)  resulting  from  the  operation  of  the  system. 
Figures  1  and  6  show  locations  of  the  stacks  and  the  15  receptors 
for  which  data  is  available.   As  previously  mentioned,  data  was 
taken  at  each  receptor  as  a  continuous  function  of  wind  azimuth 
for  at  least  one  and  usually  two  wind  speeds  (11.3  fps  and  28.3 
fps) .   The  speed  of  11.3  fps  was  chosen  because  worst  case  concen- 
trations were  found  around  this  speed  in  the  parametric  study; 
28.3  fps  is  approximately  the  annual  mean  gradient  wind  speed  for 
Boston,  and  is  thus  typical  of  what  one  would  expect  on  a  routine 
basis.   The  data  reported  reflects  the  peak  values  observed  at 
any  of  these  receptors  for  a  complete  360  degree  wind  rose.   The 
parametric  trends  shown  in  Section  2   (i.e. , Figures  3  through  5) 
would  apply  qualitatively  to  the  data  summarized  in  this  table. 

Table  1  summarizes  the  features  of  the  full  scale  ventilation 
system  and  of  the  CO  "source"  emission  rates  estimated  for  the 
worst  case.   Table  2A  summarizes  the  peak  CO  concentrations 
increments  and  the  dilution  factors  applicable  to  each  of  the  15 
receptors  for  the  highest  concentrations  observed  with  design  (full) 
volume  flow  from  Stack  I. 


*   See  Appendix  A  for  definition. 
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Figure  6.   Receptor  Locations  Beyond  Copley  Place  Site 
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TABLE  1 


PHYSICAL  DIMENSIONS  AND  FLOW  CHARACTERISTICS  FOR  RAILROAD,  AND 
PROPOSED  COPLEY  PLACE  TURNPIKE  AND  GARAGE  VENTILATION  SYSTEMS. 


Stack  I   (Turnpike  Mainline  and  Railroad) 

Exit  location  on  Stuart  Street  at  El. 

183 

8  fans 

Full  Output  Volume  Flow:  1.23  million 

cu. 

ft. 

/minute 

Approximate  Size  of  Stack  Opening:  12  0'  x 

24' 

Stack  Air  Temperature:  ambient 

CO  mass  emission  rate  for  case  in  question: 

14.06  gm/sec  (Hancock  fans  : 

running) 

17.3  gm/sec  (Hancock  fans  o 

ff) 

Corresponding  CO  Volume  flow  rate  at  standard 

conditions 

6.78  X  10   cc/min   & 

8.34  X  10   cc/min 

Stack  II:   (Turnpike  Ramps  &  Garage) 

Located  at  SE  corner  of  retail  area  a 

t  elevat 

ion  185 

3  fan  systems  -  Ramp  B  and  D  

Central  Garage 

Marriott  Garage 

Total 

.180 
.330 
.130 

000 
000 
000 

cfm 
cfm 
cfm 

.640 

000 

cfm 

Approximate  size  of  stack  opening:  32 

•  X  28' 

CO  mass  emission  rate  for  this  case: 

6.12  gm/sec  (corresponding 

volume 

flow  rate 

or  1.14  X  10^ 

cc/min 

at  STP) 
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The  CO  increments  calculated  can  be  viewed  as  valid  for  whatever 
length  of  time  the  wind  conditions  can  be  assumed  as  stationary 
in  time   and  direction;  i.e.  if  the  wind  were  steady  for  one  hour, 
the  CO  increments  could  be  added  to  background  levels  and  other 
sources  and  compared  to  a  one-hour  standard,  etc.   Since  the 
peaks  svimmarized  are  for  the  wind  direction  producing  the  worst 
case  at  each  receptor,  these  CO  concentrations  could  be  used  to 
product  a  worst-case  impact  in  deriving  either  a  one-hour,  eight 
hour  (or  other  period)  air  quality  impact  level.   It  is  anticipated 
that  the  data  supplied  will  be  used  to  perform  such  calculations 
for  the  Copley  Place  EIR/EIS. 

Table  2B  summarizes  similar  results  for  the  operation  of 
Stack  II.   It  can  be  seen  that  the  incremental  increases  in 
carbon  monoxide  level, both  near  the  site  and  far  away, due  to  the 
operations  of  Stack  II  are  small. 

For  those  wind  directions  where  impacts  are  observed  at  a 
given  receptor  for  both  Stacks  I  and  II,  the  combined  impact  of 
the  ventilation  systems  can  be  computed  by  simply  adding  the 
concentration  levels  observed  due  to  each  stack.   If  additional 
sources  are  identified  within  the  turnpike  and  garage  system, 
then  the  dilution  factors  can  be  used  to  calculate  incremental 
impacts  due  to  the  additional  source. 
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APPENDIX  A 
CALCULATION  OF  FULL  SCALE  CONCENTRATION  FROM  MODEL  SCALE  DATA 

Introduction 

When  pollutant  concentrations  or  efflux  rates  are  measured 
or  quoted  at  either  model  or  full  scale,  and  then  converted  to 
the  other  condition,  extreme  care  must  be  taken  to  ensure  that 
mass  concentration  (or  efflux)  and  volume  concentration  (or 
efflux)  are  not  inadvertantly  confused. 

Pollutant  efflux  rates  are  usually  stated  in  grams/second, 
a  mass-based  parameter,  while  air  quality  standards  are  usually 
written  in  terms  of  volume  concentration;  i.e.,  parts-per-million 
(ppm)  . 

A  detailed  discussion  of  the  general  conversions  involved 
in  fluid  modelling  may  be  found  on  p.  141  of  Synder's  EPA  Draft 
"Guideline  for  Fluid  Modelling  of  Atmospheric  Diffusion"  EPA- 
450/4-79-016  (June,  1979). 

Concentration  Ratios ; 

The  ratio  of  mass  concentration,  C,  in  full  scale  (subscript 
f)  to  that  at  model  scale  (subscript  m)  is: 
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m   ^  f'  ^  f'     ^  m 

where   U  =   gradient  wind  speed 

H 
H  =   typical  length  dimension  r^—  =  scale  factor 

"f 

M  =   mass  flow  rate  of  pollutant  (typically  in  gm/sec) 
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The  voliime  concentration  x  ratio  is  simply 
Xf    C^  /  C. 


where  p  =  the  fluid  density  of  the  full  scale  or  model  scale 
pollutants  or  tracer  gases  respectively. 

The  mass  flow  rate  M  is  simply  pQ,  where  Q  is  the  volume 
flow  rate  of  a  particular  pollutant.   In  the  case  studied,  the 
respective  gases  are  both  at  ambient  termperature  for  both  the 
model  and  full/scale  conditions;  thus, 

Xm   "   "f   Qm  l«fj 
In  the  modelling  conducted  in  this  study,  U   =  U^,  and 
H   /H^  -    1/384;  therefore  Eq.  A.  3  simply  becomes 

X.  =  X^   (6.78  X  10"^)  X  ^ 

Ttl 

where  x   =  the  measured  CH.  concentration  in  ppm,  and 

Q^  =  6.78  X  10   cc/min  of  CO  for  Stack  I  at  full  flow; 

I  5 

Q^     =  1.14  X  10   cc/min  of  CO  for  Stack  II  at 

II 

full  flow, 

Q   ranged  between  200  and  800  cc/min  in  the  model  tests,  depending 
m     -^ 

upon  the  particular  test,  but  200  cc/min  was  used  most  often  for 
the  data  quoted  earlier.   Thus,  for  example. 
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Y   =2.3xlO^Y   for  Stack  I  with  200  cc/min  of 
'^f  in 

tracer  flowing  in  the  model. 
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Dilution  Factor 

The  dilution  factor  is  a  convenient  aid  to  understanding 
the  relative  diffusion  of  pollutants  from  a  source  (such  as  a 
ventilation  stack)  to  remote  receptors.   The  dilution  factor 
(DF  )  is  defined  here  as: 


DF   = 


^measured  -  ^background 


■^stack    -  ^background 


When  the  gas  and  surrounding  air  temperatures  are  equal  at 
ambient  temperature,  in  the  model, 

^^^  Xstack  =   ^ 
Qair  ' 

where  Q  .   is  the  total  stack  volume  flor  rate. 

The  provide  a  numerical  illustration,  the  stack  flow  for 

Stack  I  was  8.2  cfm  (2.32  x  10  cc/min)  in  the  model  and  the 

-3 
tracer  flow  rate  was  800  cc/min;  thus  x  a.   ,  =  3.45  x  10 

(or  3445  ppm)  in  this  case.   If  the  measured  concentration  at 

a  given  receptor  s  100  ppm  and  the  background  concentration 

in  the  tunnel  was  5  ppm,  then  DF  =3440.5  "^    0.0276,  or  2.76% 

at  that  receptor. 
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